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Abstract 

We explore the mechanism of transverse momentum broadening of fast quarks 
propagating in nuclei, using Drell-Yan (DY) transverse momentum distributions mea- 
sured in the experiment E866 at FermiLab with beams of 800 GeV protons. Our 
theoretical analysis is based on the color dipole approach in the target rest frame, 
which has provided a successful phenomenological description of a variety of hadronic 
reactions. The present application is relevant to the regime of short coherence length 
(SCL), where the spatial extent of the fluctuations of the projectile responsible for 
the Drell-Yan reaction is short compared to the internucleon spacing. In this limit, 
momentum broadening comes from initial state interactions and is described as color 
filtering, i.e. absorption of large-size dipoles leading to diminished transverse sepa- 
ration and hence enhanced transverse momentum. The predictions we present are in 
good agreement with the E866 data. The interactions leading to the acquisition of 
transverse momentum arise from the color-dipole cross section determined previously 
from deep-inelastic scattering on proton targets. Aside from the determination of the 
color-dipole cross section, no other phenomenological input is needed to explain the ex- 
perimental results. The mean-square momentum broadening of dileptons determined 
in a recent separate analysis of the data is likewise well described by our theory. These 
results confirm within the model studied that the origin of momentum broadening in 
DY is the color dipole cross section mediating soft initial state interactions between 
the parton of the projectile that initiates the reaction and the nucleons of the nucleus, 
as provided by the color dipole description. Predictions for broadening observables at 
RHIC are presented. 
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1 Introduction 



Experimental facilities capable of studying collisions of heavy nuclei at ultra relativistic 
energies, such as the Relativestic Heavy Ion Collider (RHIC) at BNL and the future Large 
Hadron Collider (LHC) at CERN, have been motivated by the expectation that data from 
such collisions will reveal new properties of matter under extreme conditions of temperature 
and density. In recent years, interpretations of data indeed suggest that during the early 
stages of these collisions temperatures and densities are created capable of supporting the 
quark-gluon plasma (QGP), a state of matter never before produced in the laboratory and 
one of the more dramatic predictions of quantum chromodynamics (QCD), the theory of 
the strong interactions. For this reason, there is great interest in these experiments and in 
the challenging theoretical problem of determining properties of the equation of state from 
the data. These challenges derive in part from the difficulties in unfolding the equation of 
state from uncertainties in describing the reaction itself. To help resolve these difficulties, 
simpler nucleon-nucleus and deuteron-nucleus experiments, where the QGP is not produced, 
have been performed under conditions of energy and momentum transfer similar to those 
encountered in the nucleus-nucleus collisions. 

Some of the most important issues that can be resolved in nucleon-nucleus reactions are 
related to parton propagation, for example, deciding how quarks and gluons loose energy and 
acquire transverse momentum as they traverse nuclei under various conditions of coherence. 
The issues of parton propagation can be cleanly studied in Drell-Yan (DY) reactions in 
which lepton pairs produced with hadronic beams are detected under controlled kinematic 
conditions. The important characteristic controlling nuclear effects in the DY reaction is 
the coherence length, or equivalently the coherence time, which is the lifetime of the hard 
fluctuation containing the lepton pair. In the regime of short coherence length (SCL), which 
is of particular interest in this paper, the hard reaction occurs inside the nucleus over a short 
time scale, thus enabling one to identify the stage of initial state interactions , i.e. those 
preceding the hard interaction. Thus, the DY reaction in the SCL is clean in the sense that 
the produced lepton pair hardly interacts at all on its way out of the nucleus, conveying 
undisturbed information about how the parton of the projectile that initiates the reaction 
propagates in the nucleus before the production of the lepton pair occurs. The coherence 
length increases with energy, and in the regime of long coherence length (LCL) one can 
no longer separate the stages of initial and final state interactions. In this case the whole 
nucleus acts as a single scattering center. 

Notice that only in the LCL regime one can rely on QCD factorization and translate 
multiple interactions of the projectile partons in the nuclear target into a modification of the 
nuclear parton distribution function (PDF). In the SCL regime the initial state interactions 
are known [1] to break factorization. 

To be useful, it is necessary to have good statistics under the appropriate kinematical 
conditions. FermiLab data in p + A collisions [2-5] are useful in these regards and have 
been analyzed [5-8] as a source of information about both quark energy loss and transverse 
momentum broadening in nuclei. 

In this paper we are particularly interested in momentum broadening of quarks prop- 
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agating in nuclei. Theoretical frameworks that have been used to examine momentum 
broadening in various contexts include the parton model [9-11], perturbative QCD factor- 
ization [12], diagrammatic multiple scattering [13], and the color-dipole approach [14-17]. 
All these theories depend on a phenomenological characterization of the strong interaction 
to account for non perturbative physics that at present can not be calculated directly from 
QCD. Accordingly, most require a many-body transport coefficient that is fitted to nuclear 
scattering data. An exception is the color-dipole approach, where the non perturbative 
input is fixed at the two-body level by deep-inelastic scattering on protons. 

Fixing the the non perturbative input at the two-body level is advantageous when the 
theoretical formulation is able to use this information as the basic building block for pre- 
dictions of projectile-nucleus scattering. This is precisely what is done in the color-dipole 
approach, where nuclear reactions are calculated, using Glauber-Gribov theory [18, 19] and 
its generalizations [20], in terms of a color-dipole cross section fit to the deep inelastic 
data. Consequently, once confirmed under more highly controlled experimental conditions, 
the color-dipole theory can be applied with a high level of confidence in more complicated 
situations where the transport coefficients may not be able to be determined because the 
required nuclear data is unavailable. Other desirable features of the color-dipole approach 
include the fact that since it does not entail a twist expansion, all higher twist terms are 
naturally incorporated. Another is that it does not need to be corrected for next-to-leading 
order effects since gluon radiation is already included in the dipole cross section. 

The most important advantage of the color-dipole approach is that nuclear effects are 
predicted rather than fitted to data. In this work we will implement its description of mo- 
mentum broadening as developed in Ref. [15], and compare its predictions to the DY data 
mentioned above, with specific interest in refining our understanding of how partons propa- 
gate in nuclei. The data consist of ratios of differential transverse momentum distributions 
(Tijy(Pt) for a variety of nuclei [2-5], 



of the dilepton pair has been of special interest in previous studies of parton propagation. 
Extraction of the A-dependence of A (p^) was attempted in Ref. [7], using the momen- 
tum distribution ratios from the E772 experiment [2]. Unfortunately, A (p^) can not be 
uniquely determined from the ratios, requiring instead the individual cross sections appear- 
ing in Eq. (1T1). A number of puzzles have arisen from the results of this analysis, including 
an apparent inconsistency with theoretical expectations [8,21]. The situation remains un- 
resolved, but a reanalysis [8] of the E772 data along with more recent E866 data [5, 22] 
shows that the experimental errors on A (p^) were underestimated in the original analy- 
sis, suggesting that there may in fact be no contradiction between theory and experiment. 




(1) 




(2) 
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Because of this impasse, it is particularly interesting to ask whether the theory describing 
these momentum distribution ratios agrees with the phenomenological value of A (p^) that 
has been determined from them. Thus, an additional goal of our work is to determine the 
A (pr) corresponding to the momentum distributions predicted in the color-dipole approach 
and to compare them to the recent phenomenological analysis given in Ref. [8]. 

Our paper is organized as follows. In Sect. 2 we introduce the color-dipole model 
applied to the transverse momentum distributions for the DY reaction on individual hadrons. 
Section 3 discusses the extension of the theory to DY reactions on nuclei and the way in which 
the nuclear medium modifies these distributions. Section 3.1 addresses the important issue 
of coherence length, and Sects. 3.2 presents the theory and its implementation as applied to 
DY on nuclei. Section 3.3 gives numerical results for the ratios of momentum distributions 
in comparison to experiment, and Sect. 3.4 gives the corresponding results for A (p\). 
Finally, in Sect. 3.5, we present predictions for RHIC. Various technical details justifying 
the approximations made when evaluating the theory are relegated to the Appendix. Section 
4 summarizes the paper and presents our conclusions. 



2 Drell-Yan reaction on a nucleon in the target rest 
frame 

The color dipole approach in the target rest frame provides an attractive means for inter- 
preting DY data, for reasons discussed in the Introduction. The color dipole approach was 
originally proposed in [23] for hadronic interactions, and it was applied to deep-inelastic 
scattering (DIS) in [24]. Extension to Drell-Yan was developed by Kopeliovich [25] and sub- 
sequently by Brodsky et al. [26]. According to this picture, the Drell-Yan reaction occurs in 
the target rest frame when a fast projectile quark scatters off the gluonic field of the target, 
radiating a massive photon 7* of mass M as shown in Fig. [U In this frame, the 7* is a 
constituent of the projectile fluctuations, which are "frozen" by time dilation for a length 
of time t c , the coherence time, given by the uncertainty relation, 

__2E ± _ 

tc ~ M\ - ml ' K ) 

qll Q 

where E q and m q refer to the energy and mass of the projectile quark and is the square 
of the effective mass of the fluctuation 



with a the fraction of the light-cone momentum of the incident quark carried by the lepton 
pair and kj- the transverse momentum of the lepton pair. The DY reaction occurs when 
the 7* of the fluctuation is released by an interaction between one of the constituents of the 
fluctuation and a target nucleon. Subsequently the 7* decays into the observed Drell-Yan 
dilepton pair 11. 
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Figure 1: In the target rest frame, DY dilepton production looks like 
bremsstrahlung. A quark or an anti-quark from the projectile hadron scat- 
ters off the target color field (denoted by the shaded circles) and radiates 
a massive photon, which decays into the lepton pair. The photon decay is 
not shown. The photon can be radiated before or after the quark scatters. 



The DY cross section on a nucleon is described by 
d 2 a DY a 



M 1 



dM 2 dx\ 6ttx 



- I' dx q F*(x q ) [ d 2 p\^(a,p)\ 2 a qq (ap), (5) 

1 J Xl J 



where u qq {r) is the universal color-dipole cross section, p is the transverse distance between 
the 7* and quark in the fluctuation, and \1/ represents the light-cone distribution amplitude 
for the incident quark to fluctuate into a quark and the 7*. In DY, the color dipole consists 
of the quark before and after the release of the 7*, whose impact parameters differ by ap. 
In Eq. ([5]) x q (xi) is the fraction of light-cone momentum of the incoming hadron h, for us a 
proton, carried by the quark (lepton pair) and F£ is the average quark distribution function 
of the incident proton, 

=£[?/(*«) +?/(*«)]• (6) 
/ 

The light-cone momentum fraction a is given in terms of these variables as a = x\jx q . 

Deep-inelastic scattering from nucleons at HERA (DESY) has been used to fix models of 
the color-dipole cross section. The HERA data suggests that the color dipole cross section 
saturates at large tt, i.e. approaches a constant value. Golec-Biernat and Wiisthoff (GW) 
chose a simple saturated form for the color dipole cross section [27] and found that the 
following expression, 

a qq (r T} x)=a (l-e~ r ^ R2 ^) , (7) 



parametrizes well these data. For small ry, 



° qq (r T )^Cr 2 T (8) 

where the factor C = ao/B^(x) depends on Bjorken x. For us, the most important point 
is that once the color dipole cross section has been fitted to experimental DIS data on 
the nucleon, the color-dipole formalism makes predictions for Drell-Yan both in nucleon- 
nucleon [28] and nucleon-nucleus scattering with no further adjustment of parameters. 

The GW theory is known to describe DY on a nucleon, in the color dipole approach, only 
at small Bjorken x, and specific calculations [28] for X2 < 0.1 have shown it to agree well 
both in the magnitude and shape with the next-to-leading order parton model. Although 
the calculations disagreed with the E772 FermiLab data, they are in excellent agreement 
with data from the E866 experiment that has recently become available [22]. 

An expression similar to Eq. (J5j) exists for the momentum distribution. 

d 4 a DY _ oiem xi f 1 do_ Fh (xx \ da(qN -> Y x ) 
dM 2 dxpd 2 pT 3nM 2 xi + x 2 J X1 a 2 q \ aj d\nad 2 pT 

The differential cross section for heavy-photon radiation in a quark-nucleon collision was 
derived in [30], 



da(qN -> 7 *A) 



dlnad 2 k T (2tt) 2 
1 



J d 2 p x d 2 p 2 exp[i/c T • (pi - p 2 )]^*» q (a, pi)* 7 * 3 (a, p 2 ) 



x 



2 



a q q(api) + a qq (ap 2 ) - a qq {a\pi - p 2 \) . (10) 



This is described in detail in Ref . [29] . In this case, the transverse momentum pt of the DY 
pair is measured. The DY variables X\, M 2 , x 2 , and k T satisfy the equations, 

M 2 + k 2 T 

XiX 2 = , (11) 

s 

xi - x 2 = x F , (12) 

where Xp is the Feynman variable. The variables x\ and x 2 are interpreted in the parton 
model as the Bjorken variables of the annihilating quark and antiquark. 

Since the virtual photon of Eq. (TTTJT) originates in a projectile fluctuations at the hard 
scale, the quark distribution function over which this cross section is averaged to get the 
Drell-Yan reaction should be taken at the same scale. Additionally, in principle the distri- 
bution function should be of leading order since gluon radiations are included in the dipole 
cross section. In practice, the cross section ratios in which we are interested are rather sen- 
sitive to the details of the choice of the quark distribution function, and for our calculations 
we use a simple phenomenological choice [31]. 
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3 Medium effects in DY reactions in the target rest 
frame 



When the scattering takes place on a nucleus, multiple interactions with target nucleons 
can give rise to various medium effects. The coherence length determines how these arise, 
and it is useful to distinguish two limiting cases in which the theory simplifes. One is the 
SCL limit, reached when the coherence length £ c = c(t c ) becomes much smaller than the 
interparticle spacing d, l c <C d (in a heavy nucleus, d ~ 2 fin). The other is the LCL limit, 
reached when £ c 3> Ra, where Ra is the nuclear radius. 

One of the important medium effects is shadowing. In the SCL there is no shadowing 
because the duration of the fluctuation is so short that the constituents of the fluctuation 
have no time to multiply interact with the medium. In the LCL, which is applicable to 
reactions at the LHC and at RHIC under certain kinematic conditions, there is maximal 
shadowing. 

Interactions of partons in the nuclear medium also give rise to parton energy loss and 
to transverse momentum broadening, to which ratios of nuclear cross sections R A I A are 
sensitive. In DY, these effects arise in the SCL as partons from the incident hadron un- 
dergo soft interactions with nucleons of the nucleus before experiencing a hard interaction 
that immediately liberates a 7* from one of its short-lived fluctuations. For the LCL, the 
nuclear interactions that lead to shadowing are also a source of parton momentum broad- 
ening. Nuclear effects giving rise to momentum broadening and shadowing are believed to 
arise predominantly from the same color dipole cross section c^ry) that mediates the DY 
reaction on a nucleon. 

Cases intermediate between the SCL and LCL are generally more difficult to model. 
Several methods have been proposed for this. The most strict one is the Green function 
method [20]. Another one is a simple interpolation between the SCL and LCL limits using 
the square of the longitudinal form factor F\(q c ) expressed in terms of the longitudinal 
momentum transferred in the reaction, q c = l/l c . The latter method was recently applied 
to determine quark energy loss from ratios of pr-integrated nuclear DY cross sections, using 
E772 /E866 FermiLab data at E p = 800 GeV □. The kinematics of this data corresponds, 
for the most part, to the SCL regime; therefore, initial state interactions are predominantly 
soft and they decouple from the hard DY reaction. 

Predictions for nuclear broadening in DY reaction based on the theory [30] for the long 
coherence length limit were recently given in Ref. [29]. Here we concentrate on the SCL 
limit for which the theory to describe quark transverse momentum broadening has been 
presented in Ref. [15]. 

^dn this analysis energy loss was assumed to be independent of the parton energy in accordance with 
perturbative QCD calculations. However, for soft constituent partons an effective energy loss turns out to 
rise with energy [32]. Nevertheless, this fact does not affect the results of the analysis [6], but only their 
applications at higher or lower energies. 
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3.1 Coherence length vs transverse momentum 

As discussed in the previous section, there are two important limits of coherence length, 
those of SCL and of LCL, and distinct mechanisms of transverse momentum broadening 
apply in each of them. Experimental access to the LCL regime requires beams of sufficiently 
high energy to allow fluctuations of the projectile to be frozen by Lorentz time dilation so 
that their spatial extent exceeds the nuclear diameter. We do not have so far any DY 
measurements to verify the mechanism of transverse momentum broadening at play in this 
limit, although predictions for p+A, d+A, and A+A collisions at RHIC at xp = X\—x 2 > 0.5 
and at the LHC are available [29]. The experimental situation is different for the SCL 
limit. The most extensive data sets we have today are the momentum distributions of 
the E772/E866 FermiLab experiments. To better define the dividing line between the two 
regimes for the momentum distributions, we examine in this section the coherence length 
as a function of transverse momentum. 

We proceed by calculating the coherence length of the fluctuation in the incident proton 
containing the 7* of transverse momentum that eventually decays into a dilepton pair of 
momentum p T as measured in the laboratory. We proceed as in Ref [6] , averaging Eqs. (0 
over a, 



J dx q F*(x q )Jd 2 k T ^i lq {xi/x q ,k T ) 2 K DY (xi/x q ) 



llq\ 

(Q = 1 , (13) 

/ d.r„ /•;;■(,■,/) j * Ilq ( Xl /x q ,k T ) 



J-l 



where K (a) is 



K DY (n \_ M 2 (x 1 ,x 2 ,k T )(l- a) , v 

[) M 2 (x u x 2 ,k T )(l-a) + a 2 m 2 q + k 2 ' 1 ' 



with 



M (xi, x 2 , hr) = SX1X2 — k T (15) 

in accord with Eq. (TlTi) . 

The normalization integral for the light-cone distribution amplitude is ultra-violet di- 
vergent (for radiation of transverse photons), implying that the vacuum fluctuations are 
overwhelmed by infinitely heavy \qj*) fluctuations. Those heavy fluctuations, however, do 
not contribute to the DY cross section, since they are too small to be resolved and dis- 
tinguished from a bare quark. Therefore, a proper weight factors in the above averaging 
procedure must involve the dipole interaction cross section on top of the light-cone distribu- 
tion amplitude. The Fourier transform of such a modified light-cone distribution amplitude 
^J, was evaluated in Ref. [6] as 



i a 2 t 2 



^Jigi®! M = 2 z a V^C(s) e ■ k T - - , (16) 



where e is the polarization of the virtual photon and r = (1 — a)M (xi,x 2 ,k t ) + a m 2 . 
The calculation was done with the simplified dipole cross section a qq (rT,x) = C{x)r, 
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We show the coherence length for the FermiLab E772/E866 experiment at s 1 / 2 = 
38.3 GeV in Fig. [2j The values of x\ and x 2 relevant to the data for E772 are in the range 
0.07 <x 2 < 0.13 and 0.33 < x x < 0.43 [35]. The E866 data covers the range .02 < x 2 < .08, 
with (xi) = 0.46 and (M) = 4.6 GeV [36]. The values of px cover the range from to 
about 5 GeV/c. Although the coherence length has been evaluated using the value of x\ 




Figure 2: (Color online) Coherence length for FermiLab kinematics. 



corresponding to the E866 data, (xi) = 0.46, £ c is rather insensitive to its precise value. 
We see that the short-coherence length limit (£ c ~ 2fm) prevails for x 2 > 0.05 and for all 
values of px beyond 4.5GeV/c. The figure is cut off when M drops below 4 GeV. We see 
that this happens for p T > 4.5 GeV at x 2 = 0.05, but for x 2 > 0.07 M < 4 GeV throughout 
the range of px shown. 

We find the results for RHIC kinematics shown in Fig. [3j This corresponds to x\ = x 2 
(xp — 0) and s 1 / 2 = 200 GeV. We see that the coherence length is less than the internuclear 
spacing for x 2 > 0.05 and for all values of pt beyond 9 GeV/c. 

3.2 Calculating momentum distribution in the SCL limit 

Calculations of the previous section show that DY production for x 2 > 0.05 at FermiLab 
in the E772/E866 experiments with 800 GeV protons and at RHIC with s 1/2 = 200 GeV 
protons at xp = explore momentum distributions in the SCL limit. In this subsection we 
calculate these distributions using the theory of Ref. [15]. We will compare these calculations 
to the E866 data in the next subsection. We do not show comparisons with the E772 data 
because the transverse momentum distributions are essentially identical to those with the 
more recent and reliable data from E866. Later on, we show predictions for transverse 
momentum observables under kinematic conditions appropriate to RHIC. 
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Figure 3: (Color online) Coherence length for RHIC kinematics. 



It was recognized in that work that at high energies, the multiple interactions of a 
quark propagating in the nuclear medium can be eikonalized and exponentiate into a factor 
containing the color dipole cross section and the nuclear thickness function. The physics of 
momentum broadening of the quark is therefore that of color filtering, or absorption of large 
dipoles leading to diminished transverse separation with distance. Then, the probability 
distribution W q (kT, s q ) = dn q /d 2 kT that a valence quark arriving at a position (B, z) in the 
nucleus A will have acquired transverse momentum becomes [15] 

W 9 {kr, s q , B, z) = J d 2 b J d 2 b' e^" S ') e"SS e - W^.^W+^-O (17) 

where s q is the square of the quark-nucleon energy in their center-of-mass, and 

2 



N 



with (r 2 m ) = 0.79 ± 0.03 fm 2 the mean-square charge radius of the proton [37], and where 
Ta(6, z) is the nuclear thickness function, 

T A (b,z) = f dzp A (b,z), (19) 

J — oo 

with p A {b, z) the nuclear density. To obtain the transverse momentum acquired by a quark 
on the nucleus, W qA (kT, s q , z), Eq. (TTTj) must be averaged over the nuclear density Pa(B, z), 
which entails 

W A {k T) s q , z) = (W q (k T , s q , B, z)) N = jJd 2 Bj dzp A (B, z)W q (hr, s q , B, z) . (20) 
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Notice that in the regime of SCL the soft initial state interactions and the hard reaction 
(DY) cross section factorize. The universal dipole cross section, Eq. (J7|), fitted to low-x data 
for F 2 (x, Q 2 ), can be appropriately applied to the DY cross section on a proton, but not 
to the soft multiple rescatterings in nuclear matter described via Eq. (jl7l) . For evaluating 
Eq. (fTTI) we therefore use the KST color-dipole cross section, which covers soft and low-Q 2 
data [39]. For such data Bjorken x is not a proper variable, since the smallness of x does not 
mean high energy if Q 2 is small. The KST cross section has the same saturated shape as 
Eq. ([7]), but with different parameters that depend explicitly on the quark-nucleon center- 
of-mass energy. For DY with a proton beam, s q is given by s q = 2xim^E p , where E p is the 
laboratory energy of the protons and ttlm is their mass. Soft initial state interactions also 
lead to gluon shadowing, which we take into account by replacing a qq — > Rg(x, A)cr qq , where 
Rg{x,A) is the gluon shadowing function calculated [20] in the Green function approach 
from the KST model. See Ref. [15] for more discussion. We stress that our theoretical 
results do not rely on any adjustable parameters. 

We express the cross section, cr|5y(a,pT), for an incident quark to produce a DY pair on a 
nucleus A with transverse momentum as the convolution of the probability W gA (kT, s g ) 
with the momentum distribution for a quark to produce a Drell-Yan pair on a proton, 
a q ^ Y (a,k T ), 

a q D A Y (a,p T ) = J d 2 k T W qA (k T , s q )a q M^PT ~ ak T ) , (21) 

To obtain the momentum distribution for an incident proton, this must be averaged over 
the quark distribution function of the incident proton, as in Eq. ([9]), 

^dy(Pt) = (^Dy(a,PT))a , (22) 

where we define the bracket {f(a)) a to denote the integral of the function f(a) enclosed in 
it over the distribution function of a quark in the incident proton, 

(f(a)) a = J dadx q F^(x q ) f{a) 5(a - xi/x q ) (23) 

The combination of momenta pr — akx in Eq. (121]) accounts for the fact that the DY pair 
carries away only a fraction a of the transverse momentum kr of a quark in the incident 
proton. 

We denote the a-dependence of <7^ y (a:, kr) by g(a, kr), which we take from Eqs. (J9j) 
and (JTOD, 

g fa } k T ) = da(qN - 7 * X) ( M^^JI^I Y 1 ; (24) 
dlnad 2 kT \ d\nad 2 kT / ' 

evaluated using Eqs. (B.10) - (B.13) of Ref. [29]. The transverse momentum distribution 
CT Dy(^r) is of course obtained from cr^ y (a, kr) by convolution, 

°dy(Pt) = (^ P y(«,PT))a , (25) 

as in Eq. flU}. Thus, 

(t% y (<x> kr) = ^D Y (k T )g(a, k T ) , (26) 
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and we see that g(a, k T ) is closely related to the probability P(a, kx) = g(ct, k T )F^(xi/a)/a 
that a quark of the incident proton carrying light cone momentum fraction a of the detected 
dilepton pair radiates the 7* of momentum /c^ that decays into the pair. 

Since the production of the DY pair is a hard reaction, occurring for M > 4 GeV, we 
have used the GW color-dipole cross section to calculate g (a, kr) in Eq. ( 1241) . This said, we 
note that the only reason we need g(a, kx) is for the purpose of averaging over a, and we 
therefore expect the transverse momentum distributions on nuclei to be insensitive to the 
details of the color dipole model used for calculating it. 

After performing averages over the nuclear density and integrating over a, the momen- 
tum distribution of a DY pair produced on a nucleus, ^^(pt), is given by 

° P dy(Pt) = J d 2 k T {{W q (k T , s g )) N a q Ma,pr - ak T )) a , (27) 
and the final expression for the nuclear ratio becomes 

R A/p (Pt) = PP \ s / d 2 k T ((W\k T ,s q )) N a'g Y (u,PT - ak T ))« ■ (28) 

a DY\PT) J 

Using Eq. ffTTl) in Eq. (121]) and changing variables, we may express the transverse mo- 
mentum distribution of a quark as 

<(«,Pt) = jJj^y 2 J d 2 r T e^e-^U(r T )al P Y (a,p T - ak T ) (29) 

where 

d 2 b dz d 2 R PA (R-b,z)e ^> e -^ gq (r T , Sq )T A (b,z) _ ^ 

The average over the nucleus is examined in the appendix , where it is shown that to an 
excellent approximation the effect of averaging Eq. (Tl7l) over (B, z) is to replace 

T A (b,z)^(T A )/2, (31) 

which we will refer to as the average thickness function prescription. The average thickness 
function is defined by 

(T A ) = lj d 2 bT{b) 2 . (32) 

We may therefore obtain <Tj^(pt) by averaging over a as in Eq. (1221) with the average T A 
prescription, Eq. (|3~TI) . 

In this case, we find a numerically convenient form for U(tt) by expanding the term 
e~ r T' R o of Eq. (0), and thus obtaining, 



U(r T ) = Ae 
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This leads to 

f d 2 kT -* 

a q n Y (a,p T ) = S (a)a q ^ Y (a,p T ) + J ——^S 1 (a,k T ,p T )a'g Y (a,k T ) , (34) 

where 

S ( a ) = e^ T ^' 2 (35) 



S\(a,hr,pT) = 



as in Eq. flA.13[) of the appendix. The L = term has been pulled out and simplified 
by noting that the integrand for this term is strongly peaked about = due to the 
fact that (r 2 m ) is large on the GeV scale characterizing (t P q Y {j>t)- Notice that in the KST 
parametrization the parameter cr = °o( s ) is energy dependent and thus gets its dependence 
on a according to the convolution in ( J2TT) . This dependence, however, has a very little impact 

on (S )a- 

Because g(a, hp) is strongly peaked around a = 1, the average (Si(a, £>t,Pt)) is dom- 
inated by the mean value of a contributing to the DY transverse momentum distribution 
on a nucleon. We call this average ai(kx), defined as 

6ti(k T ) « ^{a 2 g(a,k T )) a - (37) 

Physically, ai(kj-) is the average light cone momentum fraction of the DY dilepton pair 
carried by the parton of the incident hadron that radiates the 7* of the dilepton pair. 
Note that this average depends not only on the momentum k^ of the 7* but also on X\. 
Alternatively, we could define «i(/cr) to be the solution of the equation 

Si(ai, k T ,PT) = {S 1 (a,k T ,p T )g(a,k T )) a . (38) 

In practice there is only a percent or two difference between a\ in the two definitions. 
Because Si (a, /ct,Pt) is peaked about k T = p T , the same mean value of q:i(/ct) characterizes, 
to a good approximation, the transverse momentum distribution both on the nucleon and 
the nucleus. Accordingly ai(fcr) is found to be essentially independent of A. We also find 
that for a given x\ and M, di(kr) depends smoothly on kr- 

We show ai(£>r) for 184 W in Fig. H]for x\ = 0.46, corresponding to E866. The quantity 
ai{krp) shown there can be fitted by the function 

ai(fcr) = «i(0) + {ak T + bk 2 T )/(l + ck 2 T ) . (39) 

The parameters in Eq. §M) are «i(0) = 0.977, a = -0.0281, b = -0.0178, and c = 
0.134. This function is relatively insensitive to £2 over the range encompassed by the data. 
However, at larger fc^, g(&, &t) spreads out sufficiently in a for ai(£/r) to acquire some 
sensitivity to the interval over which a is averaged in Eq. fl371) . Consequently, for larger x±, 
Gti(kr) levels off at a smaller value. 



e ^(T A )/2 (f {Ta)/2) L ^ ^ ^ ^ + _^ 



n 2 ^ L j 1_ t 1 

" L=l ' 6<r? m > 



— e R o 6(1 ™> , (36) 
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Figure 4: The quantity oti(fcr) of Eq. ( 138]) for 800 GeV protons incident on 
184 W. 



Because of the weak dependence of o"o on s 9 , (5o(a;)) Q in nearly independent of pr- With 
(So(a)) a and 6ji(/ct) so determined, our final result for the momentum distribution is given 
by 

f d 2 kj- -* 

° P dy(.Pt) = (v q DY( a iPT))a = (S Q (a)) a a^ Y (p T ) + J j—^Si(ai(hr) , hr, Pt)ody(M • (40) 
In the next section we evaluate Eq. (1401) and compare it to the experimental results. 

3.3 Comparison to FermiLab E866 DY momentum distributions 

The theoretical calculations of this section correspond to the transverse momentum distri- 
bution ct^y(pt) given in Eq. (j27|) . The distributions are determined once the color dipole 
cross section, the nuclear density, and the gluon shadowing corrections are specified. 

The gluon shadowing correction Rg{x, A) is given in Fig. 4 of Ref. [15] as a function of 
x for several values of A. In this paper we determine Rg{x, A) for the nuclei of interest by 
extrapolating from the values given there. The value of x at which we need to evaluate Rq is 
determined by the square of the transverse momentum k 2 that a quark acquires scattering 
off the gluon cloud of a nucleon at quark-nucleon center of mass energy squared s q and 
is given by x = 4fc 2 / s q . The average value of k 2 is fixed by the position of the peak of 
the effective gluon density in the KST model of the color-dipole cross section at s q . For 
s q = 500 GeV 2 , we find k 2 = .5 GeV 2 [15] and accordingly x = 10~ 3 . Although the 
gluon shadowing correction grows with the size of the target, our calculations show that the 
momentum distribution is rather insensitive to Rg{x, A) for the SCL kinematics considered 
in this paper. 
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The results that will be shown were evaluated using Eq. (j40p . which was obtained from 
Eq. ff27j) by a series of approximations, accurate at the few percent level, and made to 
facilitate the interpretation of the results and to streamline the numerical calculation of 
them. This calculation entails integrating the product of Si(a\, pr, kx), given in Eq. ( l36l . 
and cr^f r (fc r ) over kx using the parametrizations of di(kx) and a^yikr)- We will see in 
the next section that the transverse momentum dependence of cii(fcr) and <J p ^ Y {kT) are key 
ingredients for understanding A(p^). 

We begin our study with the examination of the E866 FermiLab fixed target experiment 
with an 800 GeV proton beam. The relevant data are nuclear distribution ratios of Eq. (CQ) 
found in Refs. [4,5,36]. The cross sections appearing in these ratios are doubly differential 
distributions binned in M 2 (or x\) and x 2 , 

t4 pA 

s IjWTr <41) 

The individual cross sections that appear in these ratios were not determined as accurately 
as the ratios and are therefore not available. The data is binned in unique x 2 intervals 
within a tight range in x\, i.e. (x\) ~ 0.46 [36]. 

Of these sets, those which have 0.1 > x 2 > 0.05 correspond to the short-coherence 
length limit, (£ c ) < 2 fm, where the momentum broadening is determined by initial state 
interactions, and it is this data set to which we compare our theory. For the E866 data 
corresponding to smaller X2, the coherence length begins to exceed the internucleon spacing 
at the center of the nucleus. The coherence length remains at the same time smaller than 
the nuclear diameter for these data sets, so these data correspond to a situation where 
the SCL and LCL mechanisms are both relevant, requiring Green's function techniques or 
interpolation, as discussed earlier. 

We would like to avoid data for which the LCL mechanism begins to play a role. The 
reason is that the mean-square transverse momentum corresponding to the LCL mechanism 
is poorly defined [29]. Although our recent analysis in Ref. [8] found no evidence for the 
mean-square momentum to depend on coherence length, the large statistical errors do not 
rule out as much as a factor of two variation among the various data sets. Thus, we restrict 
our study to those momentum distribution ratios that correspond to short coherence length 
with 0.1 > x 2 > 0.05 and for which the measured broadening of dileptons is related to soft 
initial state interactions and broadening of the projectile quark. 

For a^yikr) we have chosen a parametrization similar to that in Ref. 

1 + 4 e -4MS 

°&a*) = N \ 1+ M®.' ■ < 42 > 

where the factor in the numerator is introduced to describe a possible forward minimum in 
the cross section indicated by data [22] . It was found that the shape of the px dependence 
does not vary with dilepton mass and x 2 . A global fit to Drell-Yan data from pp collisions 
measure by the E866 experiment for different M and x 2 bins [22], with common parameters 
n, Xi and different normalization factors N, led to quite a good description [40], x 2 /d.o.f. = 
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Figure 5: Comparison of the theoretical prediction for R w ^ B£ (pt) vs. px 
(in GeV/c) with experiment for rr 2 — 0.05. Data are from the FermiLab 
E772/E866 collaboration [4, 5, 36]. 



1.5. The values of the parameters are Ai = 0.74 ± 0.12; A 2 = 0.59 ± 0.06; A 3 = 3.03 ± 0.25; 
n = 6.24 ±0.74. 

Notice that Eq. (125]) is not the same as ct^y(pt) calculated as in Eq. (TTO"]) . By using 
the phenomenological proton-proton transverse momentum distribution we compensate for 
the fact that a portion of the transverse momentum distribution arises from the density 
matrix of the initial quark in the projectile that is not yet accounted for in Eq. (flO]) . This 
explains why the measured DY momentum distribution falls off more slowly with transverse 
momentum than Eqs. ([9]) and (fTOj) evaluated as explained in the text. 

Our prediction for R w / Bs (pt) is shown in Fig. For this calculation, the gluon shad- 
owing correction was determined to be Rq = 0.9 for 184 W and Rq = 1.0 for 9 Be. One 
sees that the theory is in quite good agreement with experiment. Notice that the details 
of the Cronin effect, including the rise above 1 for pt ~ 3 in Fig. [5j are completely ex- 
plained with no adjustable parameters. For W, the value of < T4 >= 1.35, and for Be it is 
< T A >= 0.288. We use densities from Ref. [41]. 

The importance of the momentum dependence of a^kr) is illustrated in Fig. 0, where 
we compare the calculation of R w l p with a,\ — ► a,\{kr = 0) with the same quantity, but 
now with q:i(/ct)- We see that retaining the dependence on kx makes R w l p drop somewhat 
faster at larger values of pt- 

Our prediction for R Fe / Be (pT) is shown in Fig. [71 For the calculation, we found the 
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Figure 6: R w l p (p T ) for 184 W with ai = oii(p T ) (solid curve), and with 
on = oii(pt = 0) = 0.9775 (dashed curve). 



gluon shadowing correction to be Rg = 0.95 for Fe. We find comparable agreement for 
both R w ' Be and R Fe ' Be . 

1.3 \ I ■ ■ r 1 '[ ■ i 
1.2 
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Figure 7: Comparison of the theoretical prediction for R Fe l Be (j>T) vs. pr 
(in GeV/c) with experiment for %2 = 0.05. Data are from the FermiLab 
E772/E866 collaboration [4, 5, 36]. 
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3.4 Transverse momentum broadening 

The transverse momentum broadening for a dilepton produced on a nucleus is defined in 
Eq. (j2j) where the mean momentum is defined as, 



JdpWDYM 



(Pt) = r , 2 pA } x • ( 43 ) 



The dependence of the broadening A(p|.) on A for the production of heavy lepton pairs 
has been studied experimentally using the E772 data in Ref . [7] , which was revised recently 
in [8] with the inclusion of the E866 data. It was found that for large A 

&( P 2 T )i p = D(A/2)^GeV\ (44) 

where the factor D was found to be for the two sets of data, 

D(E772) = (0.029 ± 0.008 + 0.009) ( GeV/c) 2 ; 
L>(£866) = (0.059 ± 0.009 + 0.01) ( GeV/c) 2 . (45) 

The last correction here is the systematic error, with a sign known to be positive. This 
correction, if added, must be applied to both E772 and E866 simultaneously. It comes from 
the theoretical uncertainty related to the lack of data at large px- 

We next calculate A(py) in the theory just described. Because the evaluation of (p\) in 
Eq. (j43j) entails an integration over all transverse momenta, we are lead to delta functions 
that require evaluations of various quantities at Tt = 0. For this reason, we begin by 
transforming to coordinate space. From Eq. (T29]) we easily find that 



<y(Pt) = {a < S Y {a,r T )) a = j j dh T e^ T e~^U{ar T ){a < E Y {a } r T )) a (46) 
where U(tt) is given in Eq. (130!) and 

« y (a, r T )) a = a p D p Y (r T ) = J ^e~^ (a% Y (*,p T )) a • (47) 

Multiplying Eq. ( H6l) by p\ and integrating over all pr, we find, using the fact that ^(r^) 
is proportional to r\ at small r T (see Eq. (jSj)), 



(Pt) A = (Pt) N + J d " b j dz j d2 MR - b, z)e-^T A (b, z) , (48) 

where 



{ " t) - *M |r - + 3(^> (49) 

is the mean transverse momentum of DY reaction on a proton target, and where 

(a 2 ) = T~qp~f — ^("MfyK r T = 0)) a ■ (50) 
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The quantity (a 2 ) in Eq. (ISU|) is the fraction of the light cone momentum carried by 
the parton of the projectile that radiates the 7*, averaged over the transverse momentum 
for fixed X\, as can be seen by transforming back to momentum space. Using the Fourier 
transform in Eq. ( 1471) and the definition of uf^a, kx) in Eq. ( l26l) . with the quantity ai(pr) 
shown in Fig. [6] as defined in Eq. ( |37|) . Eq. ( 1501) may be equivalently written as 



$ d 2 p T (TDY(PT)Ul{PT) 



(a 2 ) = J r±»uY\r±j~L\ru ^ 

J d 2 p T a p ^ Y (p T ) 

We find that (a 2 ) ~ 0.86. This value is slightly smaller than the previous estimates in 
Ref. [29] ((a 2 ) ~ 0.9) because the integrand in Eq. (}5"Tj) peaks at slightly higher p T due to 
the fact that the phenomenological momentum distribution of the proton (t p ^y{pt) is softer 
than the theoretical distribution used in Ref. [29] (see discussion below Eq. 125"]) ). 

The first term in Eq. (j4"9]) corresponds to the contribution to the mean transverse mo- 
mentum arising from the hard interaction, and the second term is the contribution arising 
from the primordial transverse momentum of a quark in the projectile hadron. The contri- 
bution to (pt) N arising from the hard interaction is formally infinite. As long as there is no 
shadowing, as in the SCL, this gets subtracted to obtain A(p^), 

A(p 2 T ) ee (p 2 T ) A - (p 2 T ) N = AC ^ 2) J d 2 b JdzJ d 2 Rp{R - b, z)e-^T A (b, z) . (52) 

This cancellation is not complete for the LCL case [29], but we have avoided the associated 
complications by restricting our attention to 22 > 0.05, where the SCL dominates. Using 
Eq. (1521) and dropping the center-of-mass density contribution (the second term in Eq. (1A.8I) 
of the appendix), which is only about a 3% correction here, we find 

A(p 2 T ) = (a 2 )C(T A ) , (53) 

where we have used the average of T A as given in Eq. (132|) . We recall [15] in passing that 
Eq. (1531) for broadening of a dilepton pair produced in DY is closely related to the broadening 
experienced by a quark propagating through nuclear matter, differing importantly by the 
following considerations. First, on average a quark in the SCL limit passes through only 
half the nuclear thickness before it produces the dilepton pair (thus Eq. (1531) is missing the 
factor of 2 that appears in Eq. (20) of Ref. [15]). Secondly, the quark gives up only a fraction 
a of its transverse momentum to the dilepton pair, which accounts for the extra factor of 
(a 2 ) in Eq. ([53]). 

We may now compare with the experimental result in Eq. (jSj). Using values of C = 4.475 
corresponding to the KST model used for our calculations of the momentum distributions, 
and taking the uniform density model for the nuclear density, 

(T A ) = ^ Po R A , (54) 

where po ~ 0.16/m~ 3 is the central density of heavy nuclei and R A ~ l.lAs, we find from 
theory that 

A 0°T>Th ~ 0.049(A/2)^ GeV 2 , (55) 
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or D = .049. This agrees with the central value of D(E866) given in Eq. (1451) to within 
the statistical plus systematic errors and is only slightly outside the errors on D(E772). It 
was concluded in Ref. [8] that large systematic errors must be assigned to the value of D in 
the original E{772) analysis, and our results in Eq. ( |55l) are certainly consistent with this 
observation. 



3.5 Predictions for transverse momentum observables in the SCL 
at RHIC 

Since momentum broadening in DY uniquely characterizes the momentum broadening of 
the quark in the initial state, a measurement of DY at RHIC would complement information 
obtained from momentum distributions obtained from production of pions, which are addi- 
tionally sensitive to final state interactions. As seen in Fig. 3, the SCL prevails at RHIC for 
DY in proton-nucleus collisions when x^ Hlc > 0.05, for a range of transverse momentum 
up to several GeV/c. This regime at = is not easily accessed experimentally 

with current detectors at RHIC, but for the reason given before, it would be interesting to 
measure. The light-cone target-rest-frame formulation provides a theoretical framework for 
calculating the transverse momentum observables for RHIC, and we next use this theory 
for making predictions. 

In dA collisions at the RHIC collider beams of deuterons are incident on nuclei at an 
energy of 100 GeV/N, corresponding to collisions between the projectile proton with a 
nucleon in the nucleus at y/s = 200 GeV. In the target rest frame, this corresponds to 
protons incident on the nucleus with laboratory energy of about E^ HIC = 21 TeV. However, 
for Xi HIC = .05, corresponding to the smallest value of Bjorken x\ at which the SCL limit 
prevails, the quarks of the beam proton that radiate the dilepton pair have an average 
energy of only about < E^ HIC >= Xi HIC E^ HIC = 1.1 TeV, which is not so different from 
the average quark energy in the E866 FermiLab experiments. As indicated earlier, the 
FermiLab E866 experiment involves larger x\ ~ 0.46, which is mainly responsible for the 
small differences we predict. 

For RHIC, we will do our calculation for the lowest value of x\ for which the SCL prevails, 
x\ = X2 = 0.05. For this value of x\, we find that the dependence of ai(kx) is given by 
Eq. (J3SD with a(0) = 0.9513 a = 0.0168, b = -0.137, and c = 0.197, corresponding to a 
somewhat stronger fall-off of oti{hr) with kr than we found for E866. 

For gluon shadowing, the larger RHIC energies do not lead to much different Rq(x,A) 
from the values found for heavy nuclei at FermiLab. For determining the effective value 
of x = 4:k 2 /s q [15] for R G (x,A), we take s q = 2mE^ HIC = 2000 GeV 2 , corresponding to 
k 2 = .6 GeV 2 . This gives x = 1.1 10" 3 and R G = 0.91 from Fig. 4 of Ref. [15]. 

The transverse momentum distribution for protons on gold, according to our theory, 
is shown in Fig [SJ We see a clear Cronin peak at pt = 2 GeV of about the same size 
seen in DY ratios at FermiLab and in pion production at RHIC. The ratio to protons 
was calculated assuming that the transverse momentum dependence of DY production on 
protons is the same as measured in E866 at FermiLab and given in Eq. fT42|) . One should of 
course use one more appropriate for the RHIC kinematics, but this is presently unknown. 
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Using the distribution measured in E866 is reasonable since the kinematics are similar 
and presumably differences will tend to cancel in taking the ratio, but still there may be 
quantitative differences when the result is calculated with a known momentum distribution 
on the proton. 
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Figure 8: Theoretical prediction of R Au ^ p (pt) vs. pt (in GeV/c) for DY 
at RHIC, yfe = 200 GeV for x F = (x 2 = 0.05). 



The value of A(p|r) corresponding to the distribution given in Fig. [8] is given by Eq. fl53|) 
with (a 2 ) = 0.67, C(s q ) = 6.04, and (T4) = 1.5. The value of (a 2 ) is somewhat smaller 
than it was for E866 since ai(/cr) falls off more quickly in k F because of the smaller x±. 
This decrease (accidentally) almost completely compensates the increase in the transverse 
momentum acquired by the quark before it produces the dilepton pair, which is determined 
by C(sg) from the KST color-dipole cross section at = and has a weak logarithmic 
increase with energy. We find 

HpDth ~ 0.042(A/2)^ GeV 2 , (56) 

corresponding to D = 0.042, very close to our theoretical value of D corresponding to E866 
and given in Eq. (1551) . 

4 Summary and conclusions 

The color dipole approach formulated in the target rest frame provides a natural frame- 
work for examining nuclear modifications to parton propagation in collisions of hadronic 
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projectiles on nuclei. In this paper, we are interested in transverse momentum broadening 
of fast quarks propagating in nuclei, which we have examined in a particular kinematic 
regime using the color dipole approach. The color-dipole approach has predictive power 
since it relates transverse momentum observables in nuclear collisions to the color-dipole 
cross section determined from deep- inelastic scattering on a nucleon target. Once the color- 
dipole cross section is determined, predictions for these observables can be made with no 
parameters to be adjusted [15]. Transverse momentum broadening in Drell-Yan production 
has been recently measured with good accuracy in p + A collisions in the E866 experiment 
at FermiLab [4, 5, 36] with E p = 800 GeV. In this paper, we have applied our theory to the 
FermiLab data. 

In the short coherence length limit, i c <C Ra, energy loss and momentum broadening 
are mediated by multiple interactions of the incident quark with target nucleons before the 
emission of the 7*. The SCL limit is appropriate for the FermiLab data and for application 
of the theory given in Ref. [15]. The transverse momentum distributions we calculate from 
this theory, presented in Sect. 3.3, are in good agreement with the measured ones. The value 
of A(py) determined from the same theory gives the result presented in Eq. ( |55|) . which is 
comparable to the empirical determination of the same quantity [7,8]. 

We have thus confirmed within the model studied that transverse momentum broadening 
in the SCL arises from soft initial state interactions experienced by the parton initiating the 
hard reaction, as described in the color-dipole approach. In addition, our results support the 
validity of the method of determining the mean-square transverse momentum distribution 
as given in Ref. [8]. 

Predictions for transverse momentum observables at RHIC in the short coherence limit 
have been presented, complementing our earlier study [29] of broadening in the long coherent 
length regime. 

5 Acknowledgements 

We would like to acknowledge our collaborators I.K. Potashnikova, A.V. Tarasov, J. Raufeisen, 
M. Leitch, J. Moss, and P. McGaughey. This work was supported in part by U.S. Depart- 
ment of Energy, the Research Ring "Center of Subatomic Studies" (Chile), Fondecyt (Chile) 
grant 1050519, and by DFG (Germany) grant PI182/3-1. 



22 



Appendix A: Average over nuclear density 



In this appendix, we are interested in computing the average over the nuclear density in the 
expression for the momentum distribution of a DY pair produced on a nucleus, <y p DY {PT)i 

° P dy(Pt) = J d 2 k T ((W q (k T , s q )) N a q Ma,p T - ak T )) a . (A.l) 

For the purpose of evaluating the nuclear average, we note that the a dependence in 
<T£y(a:,pr), as given by Eq. (fTUj) . is sharply peaked around a = 1. This means that it 
is justified to replace the more slowly- dependent functions of a by averages. One source of 
a dependence is the quark energy s q that appears in W q {k,T,s q ) through the KST color- 
dipole cross section. In the target rest frame, s q is related to the quark energy E q by 
s q = 2m p E q , with E q = E v x\ja. As the color-dipole cross section depends only weakly on 
s q , the average is insensitive to the specific value, which on average is about E q w E p X\. 
The remaining a dependence is the explicit factor appearing in Eq. flA.lj) . 

Then, the momentum distribution o"g^(p T ) may thus be expressed, to an excellent ap- 
proximation, as 

<y(Pt) « Jd 2 k T (W«(k T ,s q )) N (a% Y (a,PT-ak T )) a (A.2) 
« Jd 2 k T W^k T ,7s q )a p D p Y (p T -ak T ) (A.3) 



where a\ may depend on kx, «.e. a\ = ai(kx) and where o"jrf y (/c:r) is defined as in Eq. (1251) . 
Using Eq. f|T7|) in Eq. flA.31) and changing variables, we may obtain expressions for u P j^y{p t ) 
and U{r T ) identical to those in Eqs. (I29|30l) . 
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<y(Pt) = jJ^J d 2 r T e rk ^e-^U(r T )a p D p Y (PT ~ ak T ) (A.4) 

d 2 b dz d 2 Rp A {R-b,z)e~^^e-^ rT > s '> )TA(b > z) . (A.5) 

except that a± and s q replace a and s q . 

In this appendix, we calculate ct^(pt) in specific models, and show that we are able 
to get quantitative agreement with the exact expression, Eq. f ]A.5j) . if we make the average 
thickness function prescription of Eq. (|3~T!) for the model. 

One model we consider is the uniform density model, 

p A (b,z) =p 6(R -r), (A.6) 

where Rq and po are fixed to preserve A and (Ta) of a realistic density for the nucleus in 
question. We find in practice that when R and po are fixed in this fashion, the uniform 
density model can be applied to calculate cr^ for nuclei as light as 9 Be. For light nuclei, 
we use a gaussian density, 

PA(b,z)=p e^ b2+z ^ R o . (A.7) 
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where Rq and po are again fixed to preserve A and (T A ) of a realistic density. 
To calculate U(tt), we begin by writing 

Pa(-R - u, z) = p A (u, z) + 5p A (R, u, z) (A.8) 
and noting that U(tt) niay be written 

U(r T ) = U (r T ) + 5U(r T ) (A.9) 

where 



and 



U (r T ) = J d 2 b J^ A dT A e-^ rT ~ s " )TA (A.10) 
5C/(r r ) = N d 2 u dz d 2 RSp A (R,u,z)e e -^ m (r T ,s q )T A ( u ,z) _ ( A11 ) 



We refer to Eq. flA.111) as the center-of-mass correction. 

In general, U{r T ) is characterized by a set of numbers Ul, 



U(r T ) = U L e- Lr T/ R «/L\ (A.12) 



L=0 



obtained by expanding e t t/ r o in cx^ty, s q ) (see Eq. (J7J)). With this expansion, we may do 
the integrals over in Eq. (IA.5I) . leading to 



a 



pA 1 \ 1 71 U L f d k T ^1^+^- pp f _ r . , , 

V(Pt) = tL l 7 1 TJ- / -27^ e ( cm> <y{Pt ~ k T ) ■ (A.13) 



In the uniform density model, Uql, corresponding to the first term in Eq. flA.9j) . may be 
worked out analytically to give 

Uql = ^[r(L + 3, c) - T(L + 3) + c 2 T(L + 1) - c 2 T(L + 1, c)] , (A.14) 

C cr 

where 

c=~a (T A ) (A.15) 
with (Ta) defined in Eq. f l32l) . For the gaussian density, C/ol becomes 

/•OO 

U QL = A/c dx(T(L+ 1) - + l,ce- x )) (A.16) 

JO 

where now 

c=<7 (T A ). (A. 17) 

The contribution 5Ul of the center-of-mass correction for both the uniform density model 
and the gaussian model have been evaluated numerically. For the average thickness function 
prescription, Ul may be read from Eq. ( 1331 . 

U L = Ae^/ 2 (^(T A )/2) L , (A.18) 
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with no center-of-mass correction. 

We next illustrate the validity of the average thickness function prescription of Eq. ([32]) 
for a heavy nucleus (A = 184) in the uniform density model and a light nucleus (A = 9) in 
the gaussian density model. In these cases, we evaluate ctdy(j>t) as given in Eq. ( 1A.13I) for 
s q = 600 GeV and a — 1, but similar results are obtained for energies relevant at RHIC and 
for the full range of a that enters into the calculation of the nuclear momentum distribution 
in Fig. [61 

For A = 184, we take (T^) = 1.35, which corresponds to a nuclear density proportional 
to the charge density for 184 W [41]. The parameters of the uniform density (UD) model are 
then R = 6.99 fm and p = 0.129 fm~ 3 . The result of calculating R^p, including the 
contributions from Uol of Eq. ( 1A.14I) and 5U l calculated numerically, is shown in Fig. [9j In 



W/p 



. 9 




Figure 9: Nuclear ratio R w ^ p (pt) in the uniform density model calculated 
using Eq. I\A.13\) . including the center-of-mass correction. 



Fig. [EH the result in Fig. M is compared to the corresponding quantity calculated using the 
average thickness function (ATF) prescription, R^tf'i specifically, the solid curve gives 



rZ p (pt) 



iOO — V ^ /p tlJVFKVT, (A19) 



Clearly, R^jf and R^tf are equal to about one percent or better. The dashed curve in 
Fig. [10 shows the contribution of the center-of-mass correction 5R^p to R^^; specifically, 
the dashed curve gives 

RuL p (pt) 

The center-of-mass correction is clearly quite small, less than a percent. 



lOO^g^ • (A.20) 
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Figure 10: Comparison of R UL f and Ratf (solid curve); relative contri- 
bution of center-of-mass correction to R^ D P in the uniform density model 
(dashed curve). See text for details. 



For A = 9, we take (Ta) = 0.288, which corresponds to a nuclear density proportional 
to the charge density for 9 Be [41]. The parameters of the gaussian density (GD) model are 
then R = 2.23 fm and po = 0.146 fm~ 3 . The result of calculating R^ D P including the 
contributions from Uql of Eq. (1A.16I) and 5Ul calculated numerically is shown in Fig. [TTJ 
We note in passing that essentially the same results for 9 Be are found when we apply the 
uniform density model to this nucleus. The solid curve in Fig. [12] compares Rqd to the 
corresponding quantity in the average thickness function prescription, Ratf, specifically, 

m RZt(pr) -R^) . (A . 21) 

R B ot(PT) 

As in the case of 184 iy ', -R^ p and Ratf are ec L ua l to about one percent or better. The 
contribution of the center-of-mass correction 5Rq D /p to -Rgd P , specifically, 

(A.22) 

Rgd (Pt) 

is shown as the dashed curve in Fig. [12j The center-of-mass correction for 9 Be is clearly 
quite small; as for 184 M^, it contributes less than a percent to the nuclear ratio. Because 
the center-of-mass correction in Eq. (lA.lip is very small even for light nuclei, it will remain 
small for all realistic choices of density, and we therefore neglect it in our final calculations. 

The validity of the average Ta prescription of Eq. (132]) is confirmed for light nuclei 
by observing that Uo(tt), Eq. ( ]A.10[) . is preserved under this substitution up to second 
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Figure 11: Nuclear ratio R Bs ^ p (pt) in the gaussian density model calcu- 
lated using Eq. AA.13\) . including the center- of- mass correction. 



order in Ta in its multiple scattering expansion. This observation explains why the uniform 
density and the gaussian density models work essentially equally as well for 9 Be. It is 
perhaps surprising that the average thickness function prescription works so well even for 
heavy nuclei, when the higher order multiple scattering terms are important. The result is 
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Figure 12: Comparison of R G p and R 
bution of center-of-mass correction to R{ 
(dashed curve). See text for details. 



'-ATF ( s °tid curve); relative contri- 
',q^ p in the gaussian density model 
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however intuitive, and it shows that on the average the quark propagates through half the 
nucleus before undergoing a Drell-Yan reaction. Based on our comparisons of the average 
thickness function approximation to numerical calculations in models, we conclude that the 
average thickness function approximation is justified for applications to data, and we will 
use it in the main body of the paper for calculating Opy. 
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